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A self-consistent projection operator method for single-particle excitations is developed. 
It describes the nonlocal correlations on the basis of a projection technique to the retarded 
Green function and the off-diagonal effective medium. The theory takes into account long- 
range intersite correlations making use of an incremental cluster expansion in the medium. A 
generalized self-consistent coherent potential is derived. It yields the momentum-dependent 
excitation spectra with high resolution. Numerical studies for the Hubbard model on a sim- 
ple cubic lattice at half filling show that the theory is applicable in a wide range of Coulomb 
interaction strength. In particular, it is found that the long-range antiferromagnetic corre- 
lations in the strong interaction regime cause shadow bands in the low-energy region and 
sub-peaks of the Mott-Hubbard bands. 
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1. Introduction 

Single-particle excitations play an important role in condensed matter physics. They de- 
termine among others basic properties of solids such as the metal-insulator (MI) transition, 
magnetism, and superconductivity. 1 Recently developed angle-resolved photoemission spec- 
troscopy allows to observe details of the excitation spectra in various materials. 2 ' 3 The exci- 
tations are usually strongly influenced by electron correlations. Therefore various approaches 
to treat the correlations have been developed. Hubbard, 4 ' 5 for example, was the first who 
proposed a theory of the MI transition on the basis of the retarded Green function. He de- 
rived lower and upper Mott-Hubbard incoherent bands caused by strong electron correlations. 
Cyrot 6 extended the theory to finite temperatures by using the functional integral method. 
Penn 7 and Liebsch 8 developed a Green function theory starting from the low-density limit 
using the t-matrix approximation. Fulde et al. 9-13 proposed methods which use the projection 
technique. The latter describes the dynamics of electrons by means of the Liouville operator 
on the basis of the retarded Green functions. 
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In the past two decades single-site theories of excitations with use of an effective medium 
have extensively been developed. Progress was made for the MI transition in infinite di- 
mensions where the self-energy of the Green function becomes independent of momentum. 1 
Several authors 15-18 determined the self-energy (i.e., an effective medium) self-consistently so 
as to be identical with the local self-energy of an impurity embedded in a medium. The theory, 
called dynamical mean- field theory (DMFT), can be traced back to a many-body coherent 
potential approximation (CPA) in disordered alloys 19 and is equivalent to the dynamical 
CPA 20,21 used in the theory of magnetism. 22 The DMFT has clarified the MI transition in 
infinite dimensions. One of its new features is that it can describe both quasiparticle states 
near the Fermi level and the Mott-Hubbard incoherent bands. 

In the DMFT as well as the dynamical CPA one usually deals with the temperature 
Green function. Therefore one needs to perform numerical analytic continuations at finite 
temperatures, which often causes numerical difficulties in particular at low temperatures. 
We recently proposed an alternative method which directly starts from the retarded Green 
function. The method, which is called the projection operator method CPA (PM-CPA), 23 
has been shown to be equivalent to the dynamical CPA and the DMFT. 24 In the PM-CPA, 
the Loiuville operator is approximated by an energy-dependent Liouvillean for an effective 
Hamiltonian with a coherent potential. The latter is determined by a self-consistent CPA 
condition. By solving an impurity problem embedded in a medium with use of the renormalized 
perturbation theory (RPT), we obtained an interpolation theory for the MI transition. 

Although one can treat the excitations from the weak to the strong Coulomb interaction 
regime with use of the single-site theories mentioned above, the single-site approximation 
(SSA) neglects the intersite correlations which often play an important role in real systems 
such as the Cu oxcides and Fe-pnictides high-temperature superconducting compounds. In 
fact, recent photoelectron spectroscopy found a pseudo gap 25 ' 26 and a kink structure 2 ' 28 in 
cuprates which can not be explained by a SSA. 

Because of the reasons mentioned above, we have recently proposed a nonlocal theory of 
excitation spectra called the self-consistent projection operator method (SCPM). 29 It is based 
on the projection technique 30,31 and the incremental cluster expansion method. 32-34 In this 
theory, all the off-diagonal self-energy matrix elements are calculated by means of an incre- 
mental cluster expansion from a diagonal effective medium T,a a (z). The latter is determined 
by the CPA condition. We call this method here and in the following the SCPM-0, The theory 
takes into account long-range intersite correlations, which are missing in the other nonlocal 
theories such as the dynamical cluster theory 35 ' 36 and the cluster DMFT. 37 ' 38 Moreover, the 
theory can describe the momentum dependent excitation spectrum with high resolution be- 
cause it is based on the retarded Green function and because the intersite correlations are 
taken into account up to infinity at each order of the incremental cluster expansion. Using the 
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SCPM-0, we investigated the excitation spectra of the Hubbard model on the simple cubic 



marginal Fermi liquid behavior 39 and a kink structure 40 ' 41 in the quasiparticle state in the 
underdoped region. The former was proposed in a phenomenological theory and the latter 
was found in the photoemission experiments in cuprates. 27,28 

The SCPM-0 makes use of a diagonal effective medium to describe the on-site correlations 
at surrounding sites. The self-consistency between the self-energy and the medium is achieved 
only for the diagonal site. Such a treatment results in a limited range of applications. In this 
paper, we extend the SCPM-0 by introducing a new off-diagonal effective medium Y<ij a (z). 
All calculated off-diagonal self-energy matrix elements are consistent with all those of the 
medium. Such a full self-consistent projection operator method (FSCPM) should allow for 
an improved quantitative description of nonlocal excitations in solids, and should play an 
important role in the phenomena with long-range charge and spin fluctuations, in particular, 
in low dimensional systems. 

In the following subsection 2.1, we briefly review the retarded Green function and pro- 
jection technique. In the §2.2, we introduce an energy-dependent Liouville operator whose 
corresponding Hamiltonian consists of the noninteracting Hamiltonian and an off-diagonal 
effective medium. We then expand the nonlocal self-energy embedded in the off-diagonal ef- 
fective medium using the incremental cluster method. In order to obtain the different terms 
of the expansion, we need to calculate cluster memory functions for the medium. In §2.3 
we obtain them by making use of the renormalized perturbation scheme. The off-diagonal 
medium is determined by a fully self-consistent condition. We present in §3 numerical results 
for momentum dependent excitation spectra for the Hubbard model on a simple cubic lattice 
at half-filling, and demonstrate that the full self-consistency extends the range of application 
from the intermediate Coulomb interaction regime to the strongly correlated regime. Calcu- 
lated spectra modify the previous results based on the SCPM-0 so as to suppress the nonlocal 
effects in the intermediate regime. In the strong Coulomb interaction regime, we find the 
shadow bands associated with strong antiferromagnetic (AF) correlations and the band split- 
ting in the Mott-Hubbard incoherent bands. In the last section we summalize the numerical 
results and discuss future problems. 

2. Full Self- Consistent Projection Operator Method 

2.1 Retarded Green function and projection technique 

We adopt the tight-binding Hubbard model consisting of the Hartree-Fock independent 
particle Hamiltonian Hq and the residual interactions with intraatomic Coulomb interaction 
parameter U as follows: 



lattice and the square lattice. 



39-41 



Especially in the two-dimensional system, we found a 




(1) 
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i,a i,j,a 

Here e a = eq — n + U{ni^ a ) is the Hartree-Fock atomic level. Note that (• • • ) denotes a thermal 
average. The quantities eo, /i, and Uj are the atomic level, the chemical potential, and the 
transfer integral between sites i and j, respectively. Furthermore a\ a {ai a ) denotes the creation 
(annihilation) operator for an electron with spin a on site i, while = a\ a a,i a is the electron 
density operator for spin a, and dn^ = — (rii a ). 

The excitation spectra for electrons are obtained from a retarded Green function defined 

by 

Gija (t - = -iO(t - t')([a la (t),al(t')} + } . (3) 

Here 9(t) is the step function, ai a (t) is the Heisenberg representation of defined by 
o>icr(t) = exp(iHt)aier exp(—iHt) . Furthermore, [ , ]+ denotes the anti-commutator between 
the operators. 

The Fourier representation of the retarded Green function is written as follows. 1 

Gija(z) = (a\ a | j^j^a]a) ■ (4) 

Here the inner product between the operators A and B is defined by (A\B) = ([A + ,B] + ). 
Also z = to + i5 with 5 being an infinitesimal positive number, and L is a Liouville operator 
defined by LA = [H, A] for an operator A. Note that [ , ] is the commutator between the 
operators. 

Using the projection technique, we obtain the Dyson equation for the retarded Green 
function as 

G lja (z) = [(z -H - A(z)r V . (5) 
Here the matrices are defined by (Ho)ij a = e a Sy + Uj, and 

A ija (z) = (A(z)) ija = U 2 G i3a (z) . (6) 
The reduced memory function Gij a (z) is defined by 



The operator A\ a is given by 

4a = Ajrii-a , (8) 

and the Liouville operator L acting on A^ a is defined by L = QLQ, Q being Q = 1 — P. The 
projection operator P projects onto the original operator space {(aJo-)}, 

^ = El4r)(4r|- (9) 
jo- 
in a crystalline system, the Green function Gij a {z) is obtained from its momentum rep- 
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resentation Gka{z) as 

Gija(z) = ^ E G ^ Z ) eX P( ifc • " R j)) ■ ( 10 ) 

k 

Here is a number of site. The momentum-dependent Green function is given by 

Gka(z) = 1 . (11) 

z- e ka - k ka {z) 

Here e ka = e CT + £fe is the Hartree-Fock one-electron energy eigenvalue, is the Fourier 
transform of Uj , and 

Afc CT (z) = A jOa(z) exp(ik -Rj) . (12) 

j 

It should be noted that the momentum-dependent excitation spectrum is obtained from 

Pka{u) = — ImGto(z) . (13) 

7T 

The local density of states (DOS) is given by 

Pia(w) = --ImG iiff (z) . (14) 

7T 

which is identical with the average DOS per atom, p a (u>) = 1/N J2k Pkai^) when all sites are 
equivalent to each other. 

2.2 Incremental cluster expansion in the off-diagonal effective medium 

In the full self-consistent projection method, we introduce an energy-dependent Liouville 
operator L(z). Its corresponding Hamiltonian is that of an off-diagonal effective medium 
t, ija (z). It is 

L(z)A=[H (z),A] , (15) 

for arbitrary operator A and 

H (z) =H + Y^ %,*{*) a Ui* ■ (16) 

ijcr 

Defining an interaction Liouville operator L\(z) such that 

L l {z)A=[H l {z),A] , (17) 

with 

Hi(z) = U^2 Sn^Sn u - ^ E ija (z)al a a ja , (18) 

i ijcr 

we can rewrite the original Liouville operator L as follows. 

L = L(z) + Li(z) . (19) 

It should be noted that the interaction H\(z) contains the off-diagonal components 
f±ij a (z) (i / j) in addition to the diagonal ones flu a (z). Accordingly, we divide here the 
interaction Liouvillean L\{z) into single-site terms and pair-site ones. Furthermore we intro- 
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(20) 

(21) 
(22) 



After substituting the Liouvillean (19) into Eq. (7), we can expand the resolvent [z — V) 1 
with respect to L\{z) as follows. 

(z-L)- 1 = G + G TG , (23) 



T — Lj + LjGqLj + LjGqLiGqLi + 



(24) 



Here Gq = (z — Lq(z)) 1 and Lq{z) = QL(z)Q, while Lj(z) = QL\{z)Q. 

The T matrix operator may be expanded with respect to different sites as 

T = ^2 u i T i + ^ V; '- v 'i T n + ^ v i v j v kTijk H • (25) 

The single-site Tj, two-site Tjj, and three-site Tij k matrix scattering operators are obtained 
by setting the indices as (v { = 1, = 0), (i/j = Vj = 1,^^^ = 0), (i/j = vj = u k = 

!> u i{^i,j,k) = °)> and so on - Xt is 

T< = T (i) , (26) 



(27) 



^ijfc — - Ty - T jk - T ki - Ti - Tj - T k . (28) 

The operator T^ (c = i, ij, ijk, ■ ■ ■) at the right-hand-side (r.h.s.) of the above equations is 
the T matrix operators for the cluster c, i.e., 

T^ =Lf\l-G L^)- 1 . 



(29) 



Here L^f 1 = QL\ c> (z)Q and L\ c> (z) is the interaction Liouvillean for a cluster c: 

lJt ) {z) = Y,C?(z)+ ^\z). (30) 

iec (i,j)ec 

Note that the sums in the above equation are taken over the sites or pairs belonging to the 
cluster c. 

Substituting Eq. (25) into Eq. (23) we have 

G ija (z, fa}) = [a\ c I (Go + "iGoTiGo + £ W m G T lm G + ■■■ )A) a ) . (31) 

l (l,m) 
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In the incremental method, 32-34 we first consider the self-energy contribution due to intra- 
atomic excitations, i.e., 

Gg(z) = G iiff (z,vi = 1,^) = 0) = (Al \(Go + G T^G )Al) . (32) 

Next we consider the scattering contribution due to a two-site increment in Eq. (31), 

G^(z) = G ii<T (z, vi = n = 1, v mW ) = 0) = G\Z(z) + {A\ a \G m + T u )G \Al) . (33) 

This defines two-site increment to the diagonal matrix element as 

A^(z)=G^(z)-Gt(z) . (34) 

In the same way, we consider 

Gt m \z) = G$(z) + A0f!(z) + AG^\z) + (A\ c \G,{T lm + T llm )G \Al) . (35) 

Then, we define the increment for a three-site contribution. 

AGg m \z) = G^(z)-AGf!(z)-AGt?(z)-GV(z) . (36) 

The memory function Gu a (z) in Eq. (31) is then expanded as follows. 

G iia (z) = G ( tl{z) + Y J ^Gf a {z)+ ±£ E AG^ ) (z) + -.. . (37) 

It should be noted that G ii<T {z) (c = i,ij, • • • ) in Eqs. (32), (33), and (35) is obtained from 
the T-matrix T (c ^ given by Eq. (29) as 

G\%(z) = (Al \(z - L {c \z))- 1 Al) . (38) 

Here the cluster Liouvillean L (z) is defined by 

L (c) (z) = L (z)+l[ c) (z) . (39) 

In the same way, the off-diagonal memory function is obtained as follows. 

G ija (z) = GfM + y: ^(z) + \ y: e ^ m \ z) + . . . , (40 



2 



with 



AG%?(z) = G%?(z)-G ( $(z), (41) 

AGgJ m) (z) = G^ (z) - AGiff (z)- AG% m) (z)-G^(z). (42) 

When we take into account all the terms on the r.h.s. of Eqs. (37) and (40), the memory 
function Gij a (z) does not depend on the effective medium Tiij a (z). However, it is not possible 
in general to calculate the terms up to higher orders; we have to truncate the incremental 
expansion at a certain stage. In that case the memory function Gij <T (z) depends on the medium 
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^ija(z)- We determine the latter from the following self-consistent equation, 

% ja (z) = A ija {z) . (43) 

Note that the off-diagonal effective medium Xj^z) is the self-energy for the energy-dependent 
Liouvillean L(z), i.e., t ija (z) = U 2 (A\ a \ (z - Lo(z))- 1 A) c ) = U 2 (A\ C \G A] a ). Thus the 
self-consistent equation (43) is equivalent to the condition that the T-matrix describing the 
scattering from the medium vanishes according to Eq. (23): 

(Al \G TG Al) = . (44) 

This is a generalization of the CPA. 

The present theory reduces to the previous version of the nonlocal excitations (SCPM- 
0) 29 when the off-diagonal media T*ij a (z) (i ^ j) are omitted and only the self-consistency of 
the diagonal part is taken into account in Eq. (43). When in the SCPM-0 only the diagonal 
self-energy Aa a (z) is taken into account (i.e., when we make the SSA), the result reduces to 
the PM-CPA which we previously proposed. 23 



2.3 Renormalized perturbation scheme to cluster memory functions 

The incremental cluster expansion scheme given in the last subsection can be performed 
when the cluster memory function G^^z) defined by Eq. (38) is known. We obtain here an 
explicit expression for it. The Hamiltonian H^ c \z) to the cluster Liouvillean L,( c \z) in (38) 
is given by 

=H + J2 ZijoWahaja - mi VijaWahaja + U J2 Sn H Sn H ■ ( 45 ) 

ijcr ij&c cr igc 

Introducing parameters Xij a (0 < Xij a < 1), we can divide the Hamiltonian as 

H V (z) = (z) + (z) , (46) 

(z) = H + Y, ^ija{z)a\ c a ja \^{ z ) a l a j^ > ( 47 ) 

ijcr ij&c cr 

h[ c \z) = U^2 <W<W - ^2 ^2 ^ijcr%jcr(z)a] a a jcT . (48) 

iGc ij&c cr 

Here Xij a = 1 — \ij a - Note that parameters Xij a control the partition ratio of the medium 
potentials in the cluster between the noninteracting part H^ c \z) and the interacting 
part H^\z) (see Fig. 1). The Hamiltonian H^ c \z) denotes a system with a uniform effective 
medium T,ij a (z) when Xij a = 1 (i,j G c). When Xij a = (i,j G c), H^ c \z) denotes a refer- 
ence system with a cluster cavity in the effective medium and h[ c \z) denotes the Coulomb 
interactions on the cluster sites. 

According to the partition of the Hamiltonian (46), we introduce a noninteracting Liou- 
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Fig. 1. Partition of the medium on the cluster sites. The cluster Hamiltonian H^ c \z) with on-site 
Coulomb interaction embedded in the off-diagonal medium X (the left-hand-side) can be split into 
an effective Hamiltonian H^ c \z) with cluster potential AS and the remaining interaction H^ c \z). 
The parameter A can control the magnitude of the potential E on the cluster sites. Note that each 
site in the cluster may be away from the other sites of the cluster. 



villean L^\z) = QL^\z)Q and an interacting Liouvillean L,f\z) = QL^f\z)Q as follows. 

4\z)A=[H^(z),A] , (49) 



L ( f\z)A=[Hf\z),A] . 

-(c) 



(c), 



Then we have L {c \z) = Lq\z)+L\ c \z), and the cluster memory function G i C j a (z) is expressed 



(50) 



7=y( c ) 



as 



G\%(z) = (Al \(z - l£\z) - L^iz))- 1 A) ) . 



(c), 



-(c), 



(51) 



T(c), 



Note that Lj (z) in the above expression has been redefined by Eq. (50), i.e., Eq. (30) in 
which T,ij a (z) has been replaced by Xij a T,ij a (z). 

The interaction Liouvillean L : (z) expands the operator space to 

L[ c \z)\Al) = [U(l - 2<n,_ CT )) - \ lla t lla {z)] \A\ a ) - £ ^{bJJ + <^Kj) ' ( 52 ) 

(»j')ec 



\Bj ja ) = ^lja^lja(z)\a] a 8ni_ a ) 



+ h3-<T^l3-<T{z)\a\j{a]_ a ai^ a )) + A /j _ CT S /j _ (T (z)|a z T (r (5(aJ_ (T a i _ (T )) . (53) 



t S(J 



Therefore we introduce the projection operator 



(54) 



where \ic = (^i-o-)(l — (ni-a))- It projects onto the space {l^l^)}, while Q = 1 — P eliminates 

the space {l^.)}- Making use of these operators one can divide the interaction Liouvillean 
( c \ 

L : (z) into two parts. 



Lf \z) = P if \z)P + L%{z) , 



(55) 
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L%{z) = QLf\z)P + Lf{z)Q . (56) 

The first term at the r.h.s. of Eq. (55) acts in the subspace {|A[ CT )}, while the second term 
expands the operator space beyond {l^jo-)}- It should be noted that the first term at the r.h.s. 
of Eq. (56) does not vanish in general when we take into account the off-diagonal character 
of the effective medium. This point differs from the case of the SSA. 23 

Making use of Eq. (52), we obtain the expression of P lff\z)P in Eq. (55) as 



PLf(z)P = £ Kt)(4 C) W (A] a/ \ , (57) 

ijaa 



- ja ^ (l-Sa). (58) 

XiaXja' 

Substituting Eq. (55) into Eq. (51), we obtain 

G$,{z) = • (1 - 4 C) • G$ V 1 ] . . . (59) 

Here the screened memory function G 0i j o .(z) = (G )i a ja is defined by 

G~tU z ) = (4,1 {* - M*) - L%(z))~ l A) a ) . (60) 
The expression of the cluster memory function (59) is exact. The simplest approximation 

_( c ) 

is to neglect the interaction Liouville operator Ljq in the screened cluster memory function 
G 0i j a , which we called the zeroth renormalized perturbation theory (RPT-0). The approx- 
imation yields in the weak interaction limit the correct result of second-order perturbation 
theory, as well as the exact atomic limit. 

An explicit expression of the screened cluster memory function in the RPT-0 can be 
obtained approximately as shown in Appendix. 

3(c) (z)=A f dede'de"p^(X, e)p^_ a (X, e^p^jX, e") X (e, e>, e") 
° iJ %3 ° J z-e- X a t a (e, z)-e>- A_ CT £_ CT (e', z) + e" + A_ CT S_ CT (e", z) ' 

X (e, e', e") = (1 - /(e))(l - f(e'))f(e") + f(e)f(e')(l - f(e")) . (62) 

Here < X a < 1, and f(uj) is the Fermi distribution function. 

The prefactor Aij a in Eq. (61) has been introduced to recover the exact second moment 
of the spectrum. 

^ = <n_> c U-V CT > c )^ + 1 " ■ (63) 
Here the electron number for a cavity state (nj CT ) c is defined by 

(n ia ) c = J dep%i(X,e)f(e) . (64) 
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The densities of states for the cavity state are given by 

P&V) = --Im [(F^zy'+K^iz))- 1 }^ , (65) 



7T 



(F C (A, z)) l]a = F ija (z) = £ mi ^ 3 l . (66) 

— ~ (c) 

Furthermore, A CT = 1 — X a , (S {z))ij a = T>ij a (z) G c), and Y>k a {z) is the Fourier 

transform of T lija (z); T, ka (z) = £V SjOa(^) exp(zfc • Rj). 

The simplified self-energy Y, a (e,z) in Eq. (61) is calculated from Y,ij a (z) as 

S CT ( g ,z)=^ 5l e " e i ^^' ( 67 ) 

i P l e fc crJ 

where Pij(e) is the DOS of the noninteracting system defined by Pij(e) = Sfc(^)^( e — e k)(k\j) 

and p°(e) is defined by p°(e) = 1/iV <^(e - e k ). 

Finally the simplified cluster memory function in the RPT-0 is given by Eq. (59), Eq. 

(58) in which {Xij a } has been replaced by {A CT }, and Eq. (61). In the present scheme, we first 

assume T<ij a (z). Then we calculate the coherent Green function Fij a (z) (Eq. (66)), the screened 

( c \ 

cluster memory function G 0i j a (z) according to Eq. (61), as well as the atomic frequency 
matrix (Lj^)^.^, given by Eq. (58). Using these matrices, we calculate the cluster memory 
function G\f a (z) (Eq. (59)). Note that the static quantities (Aj a ,\Bj la ,) in Eq. (58) have to be 
calculated separately, for example, by means of the local ansatz wavefunction method. 1 After 
having calculated cluster memory functions, we can obtain the memory functions (37) and 
(40) according to the incremental scheme. Then we calculate the diagonal and off-diagonal 
self-energies (6). If the self-consistent condition (43) is not satisfied, we repeat the above- 
mentioned procedure renewing the medium Sj J(T (z) until the self-consistency (43) is achieved. 
When we obtain the self-consistent solution £j JO -(z), we can calculate the excitation spectrum 
Pka(u) from the Green function (11) and the DOS p a {oo) from Eq. (14). 

3. Nonlocal Excitations on a Simple-Cubic Lattice 

We present here the numerical results of excitation spectra of the Hubbard model on a 
simple cubic lattice at half-filling in the paramagnetic state in order to examine the nonlocal 
correlations in the FSCPM. As in our previous calculations, we choose the parameters A CT = 
in Eq. (61); we start from the cavity cluster state (Ay = 0) for the calculation of the memory 
function. The form (61) with A CT = reduces to the iterative pertubation scheme 43 at half- 
filling in infinite dimensions. Note that we need not to calculate S CT (e, z) defined by Eq. (67) 
as well as (AJjbJjJ in Eq. (58) when A CT = 0. 

We adopt the nearest-neighbor transfer integral t here and in the following, choose 
the energy unit as \t\ = 1. The Fourier transform of the transfer integrals is given by 
e = — 2(cos k x + cos k y + cos k z ) in the unit of lattice constant a. Furthermore, in Eqs. (37) and 
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Energy 



Fig. 2. Diagonal self-energies for U = 10 in the FSCPM (thick curves), the SCPM-0 (middle-size 
curves), and the SSA (thin curves). The real (imaginary) part in each method is drawn by the 
solid (dotted) curve. 



(40) we take into account single-site and pair-site terms only, but the latters up to the 10th 
nearest neighbors. In the FSCPM, we have to perform a numerical k integration to obtain 
the coherent Green function Fij a (see Eq. (66)). We have adopted a 80 x 80 x 80 mesh in the 
first Brillouin zone for such integration. 

In the numerical calculations of the screened memory function (61), we applied the method 
of Laplace transformation. 44 This reduces the 3-fold integrals with respect to energy into the 
one-fold integral with respect to time. 

Starting from an initial set of values {T>oj(z) = 0}(j = 0, 1, • • • , M), we calculated {Fij(z)} 
(Eq. (66)), {p^{z)} (Eq. (65)), {G^(z)} (Eq. (61)), {G\f (z)} (Eq. (59)), and {G^z)} (Eqs. 
(37) and (40)), and finally obtained Aij a (z) = U 2 Gij(z). We have repeated the same procedure 
until self-consistency (43) was achieved. 

We show in Fig. 2 the self-consistent diagonal self-energy at U = 10 and compare the 
results with the SSA and SCPM-0 calculations. The SCPM-0 shifts the spectral weight of the 
SSA towards the higher energy region. The FSCPM suppresses the amplitude of the SCPM-0 
self-energy in the high energy region. In the low energy region, we find that the imaginary 
part is close to that of the SSA, while the real part is in-between the SSA and the SCPM-0. 

The off-diagonal self-energies are shown in Fig. 3. We find again that self-energies are 
suppressed by the full self-consistency. Furthermore the self-energies rapidly damp with in- 
creasing interatomic distance. The fourth-nearest neighbor contribution and the contribution 
from more distant pairs can be neglected when U = 10, though we took into account the 
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Fig. 3. Off-diagonal self-energy A(n, z) for the n-th nearest neighbor pairs (n = 1, 2, 3, 4). 



off-diagonal ones up to 10th nearest neighbors. 

The momentum-dependent self-energies Af ca (z) are also suppressed by the full self- 
consistency. As shown in Fig. 4, in the low energy region the imaginary parts of A] C(7 (z) 
hardly depend on momentum k and are close to those of the SSA. In the incoherent region 
\u\ > 2, both real and imaginary parts show considerable k dependence. For example, the 
imaginary part of Af ca (z) shows at the T point a minimum at to 6 and a second one at 
uj rs — 5, while it shows a minimum at u = — 6 and the second one at to = 5 at the R point. 

Calculated momentum dependent excitation spectra are shown in Figs. 5, 6, and 7. For 
rather weak Coulomb interaction (U = 6), we find that the quasiparticle band is reduced by 
about 30% in width as compared with that of the noninteracting band. As seen in Fig. 5, the 
Mott-Hubbard incoherent bands appear at \u\ « 7.5 around the V and R point. The spectral 
weight of the Mott-Hubbard bands is reduced by 30% when it is compared with that of the 
SCPM-0. 

In the intermediate Coulomb interaction regime (U = 10), the quasiparticle band width 
becomes narrower and is reduced by 25% as compared with that of the SCPM-0. The spectral 
weight of the quasiparticle band around the T and R point moves farther to the lower and 
upper Mott-Hubbard bands. The latters are more localized in the vicinity of the T and R 
point than in the SCPM-0, and their peaks are reduced by 25% as compared with those in 
the SCPM-0. 

In the FSCPM scheme, we can obtain the momentum-dependent spectra in the strong 
Coulomb interaction regime (U = 16). This is shown in Fig. 7. In this regime, the quasiparti- 
cle band becomes narrower and the spectral weight becomes smaller. New excitations which 
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Fig. 4. Momentum dependent self-energies at various k points along the high symmetry line T-X-M- 
R. The real (imaginary) parts are shown by solid (dotted) curves. Note that the real part at the T 
point shows a maximum at energy u> 8 and a second maximum at uj w —3, while the imaginary 
part shows a minimum at lu « 6 and a second minimum at lu w —5. The results of the SSA are 
also shown by the thick solid and dotted curves. 

appear in this regime are the subbands at \uj\ 3.3 showing weak momentum dependence. 
Because of the strong AF intersite correlations at half-filling, these excitations may be inter- 
preted as "shadow" bands due to AF correlations, whose excitation spectra in a simple spin 



density wave (SDW) picture 45 may be given by e± DW (k) = ±ye| + A 2 . Here e& is a quasi- 
particle band, A is an exchange splitting defined by A = U c ff\m\/2, \m\ being a temporal 
amplitude of the magnetic moments with long-range AF correlations, and U c g is an effective 
Coulomb interaction. 

As seen in Fig. 7 the Mott-Hubbard bands in the strong Coulomb interaction regime 
show a weak momentum dependence around \uj\ ps ±11. It should be noted that the splitting 
between the upper and lower Hubbard bands is about 22, which is larger than U = 16 (a 
value yielding the atomic limit). In the strong Coulomb interaction limit at half-filling, we 
expect strong AF intersite correlations due to the super-exchange interaction J = 4|i| 2 /[7. 
The energy to remove (add) electrons is then expected to be eo — ^nn|<^| (eo + U + znn| «•?!), 
znn being the number of nearest neighbors (znn = 6 in the present case). Thus the splitting is 
expected to be U +2znn\J\ instead of U. This formula yields the splitting 19 instead of U = 16. 
The former seems to be consistent with the value 22. We also find additional excitations at 
\u>\ « 8. These sub-bands may be interpreted as local excitations of the lower and upper 
Hubbard bands without AF correlations. 
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Fig. 5. Single-particle excitation spectra along the high-symmetry lines at U = 6. The Fermi level is 
indicated by a bold line. 




Fig. 6. Single-particle excitation spectra at U = 10. 



The total densities of states are presented in Fig. 8 for various values of U. For weak 
Coulomb interactions U < 6, the deviation of the spectra from those of the SCPM-0 is 
negligible. When U > 6, Mott-Hubbard bands appear and the DOS deviate from the SCPM- 
0. An example is shown in Fig. 9 for U = 10. We find there that the full self-consistency 
suppresses the weight of the Mott-Hubbard bands and enhances the quasiparticle peaks. 
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Fig. 8. The total DOS for U = 4, 8, 12, 16. The DOS for noninteracting system is shown by dotted 
curve. 



Resulting DOS is in-between the SSA and the SCPM-0. When U > 10, a shadow band 
develops around \u\ = 3 as seen in Fig. 8. Furthermore, for U > 12 we find two peaks in each 
Mott-Hubbard band; one is due to the excitations without intersite AF correlations, another 
is due to the excitations with strong AF correlations. 

A remarkable point of the nonlocal excitation spectra is that the quasiparticle peak at the 



16/24 



J. Phys. Soc. Jpn. 



Full Paper 



0.16 




0.14 




0.12 




0.1 




0.08 




0.06 




0.04 




0.02 










-20 



11=10 



-15 -10 



FSCPM 

SCPM-0 

SSA 




15 



20 



Fig. 9. The DOS in the FSCPM (solid curve), the SCPM-0 (dashed curve), and the SSA (dotted 
curve) at U = 10. 



Fermi level reduces with increasing Coulomb interaction U. As shown in Fig 10, the ratio of 
p(0)/psSA(0) to the SSA, monotonically decreases with increasing U, and reaches 0.59 when 
U = 20. 

Momentum-dependent effective masses along the high symmetry line are presented in 
Fig. 11. In the FSCPM calculations, we obtained the self-consistent m& up to U = 20. We 
find that are approximately equal to those in the SCPM-0 for U < 6. But, for U > 6, 
the FSCPM enhances the mass as compared with those of the SCPM-0. Momentum 
dependence of m& becomes larger with increasing U. The mass has a minimum (e.g., 9.4 for 
U = 20) at T(0,0, 0) and R (ir,ir, ir) points, and has a maximum (e.g., 10.0 for U = 20) at 
(^vr/2,0). 

The average quasiparticle weight Z = (N^ 1 Y^k m k)~ 1 vs - U curve is shown in Fig. 10. 
Quasiparticle weight in the SSA monotonically decreases and vanishes for U c = 16.0. The 
curve in the SCPM-0 deviates upwards from the SSA. The curve in the FSCPM deviates 
downwards from the SCPM-0 beyond U « 5, and is between the SSA and the SCPM-0. The 
present result suggests that the critical Coulomb interaction for the divergence of the effective 
mass is U c (m* = oo) > 30. We want to mention that for the Gutzwiller wave function a critical 
Coulomb interaction U c (m* = oo) exists only in infinite dimensions (i.e., in the SSA). 46-48 

We have also calculated the momentum distribution as shown in Fig. 12. The distributions 
show a jump at the Fermi surface, and extend beyond the surface. The basic behavior of the 
distribution is similar to that in the SCPM-0. For U < 6, the curves (nk) agree with those of 
the SCPM-0. When U > 6, the FSCPM reduces (n k ) of the SCPM-0 below the Fermi level, and 
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Fig. 10. The average quasiparticle weight vs. Coulomb interaction curves in the FSCPM (solid 
curves), the SCPM-0 (dashed curve), and the SSA (thin solid curve). The DOS divided by the 
SSA one at the Fermi level is shown by dotted curve. Note that pssa(0) equals the noninteracting 
DOS p°(0). p(0)/psSA(0) beyond J7 C (SSA) = 16.0 stands for p(0)/p°(0). 



enhances it above the Fermi level suggesting increased localization due to full self-consistency. 



4. Summary and Discussions 

We have developed a fully self-consistent projection operator method (FSCPM) for non- 
local excitations by using the projection technique for the retarded Green function and the 
effective medium. The method makes use of an energy-dependent Liouville operator L(z) 
with a corresponding Hamiltonian for an off-diagonal effective medium Sjjo-(z). It allows us 
to calculate the nonlocal self-energy Aij a (z) by using the incremental cluster expansion from 
the off-diagonal medium. Each term of the expansion is calculated from the memory function 
for clusters in a nonlocal effective medium. The latters are obtained by the renormalized per- 
turbation theory within the RPT-0. The off-diagonal effective medium Sjj CT (z) is determined 
from a fully self-consistent condition {Sj JO -(z) = Ay cr (z)}. It is a generalization of the CPA 
as given in Eq. (44), i.e., (A^ \GqTGqA^^) = where T denotes a scattering T-matrix from 
the off-diagonal medium and Go is a resolvent for the Liouville operator Lq corresponding 
to the medium. In this way we can take into account the long-range intersite correlations as 
extensive as we want in each order of expansion. We obtain momentum-dependent spectra of 
high resolution from the weak to the strong Coulomb interaction regime. 

The present theory reduces to the PM-CPA when we omit all off-diagonal matrix elements 



18/24 



J. Phys. Soc. Jpn. 



Full Paper 



2 





1 


- 


' 11=20^"" 

18 




16 




14 




12 




10 












4 


1 


U=2 

i 



M X r MR 

Fig. 11. The momentum-dependent effective mass curves along the high-symmetry lines for various 
Coulomb interaction U. The curves in the SCPM-0 are also shown by the dotted curves for 
U = 2, 4, 6, 8, 10, and 12 from the bottom to the top. 



of the effective medium Tlij a (z)(i / j) and all the off-diagonal self-energy contributions. The 
theory reduces to a self-consistent theory (SCPM-0) when we omit the off-diagonal medium 
T,ij a (z)(i j), but take into account all the off-diagonal self-energy contribution Ajj cr (z). 

We have performed numerical calculations of the excitation spectra for the half-filled 
Hubbard model on a simple cubic lattice by using the FSCPM within the two-site approxima- 
tion. We have obtained the self-consistent nonlocal self-energy up to the Coulomb interaction 
U pa 20. We found that the FSCPM suppresses the amplitudes of local and nonlocal self- 
energy Kij a (z) for an intermediate strength of Coulomb interactions, reduces the weight of 
the Mott-Hubbard bands, and enhances the quasiparticle peaks at the Fermi level in the aver- 
age DOS when they are compared with those in the SCPM-0. Moreover the FSCPM enhances 
the momentum-dependent effective mass as compared with the SCPM-0. Thus the Z — U 
curve is located between the SSA and the SCPM-0. These results indicate that the full self- 
consistency tends to suppress the nonlocal effects found in the SCPM-0. We suggest a critical 
Coulomb interaction U c (m* = oo) > 30 for the simple cubic lattice, which is much larger 
than the SSA value U c = 16. In order to obtain the explicit value of U c , we have to take into 
account larger clusters embedded in the off-diagonal medium. 

The FSCPM enables us to investigate the nonlocal excitation spectra in the strongly 
correlated region. There the quasiparticle bands becomes narrower and their weight becomes 
smaller. We found shadow band excitations at \uj\ « 3 due to strong AF correlations and 
Mott-Hubbard sub-bands at \u\ « U/2 without intersite correlations. Moreover, we found 
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that strong AF correlations can shift the Mott-Hubbard bands towards higher energies. 

In the present calculations, we investigated the effects of nonlocal correlations assuming the 
paramagnetic state. One has to extend the calculations to the AF case in the next step because 
the ground state of the three dimensional Hubbard model is believed to be antiferromagnetic 
in general. Furthermore, the numerical results of calculations in the strongly correlated region 
presented here should be extended by taking into account the contributions from larger clusters 
in the incremental cluster expansion. Improvements of the self-energies for the clusters in the 
off-diagonal effective medium are also left for future investigations towards more quantitative 
calculations of the nonlocal excitations. 
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Appendix: Derivation of Approximate Expression for the Screened Cluster 
Memory Function (61) 

An approximate expression (61) of the screened cluster memory function in the RPT-0 
can be obtained as follows by using the momentum representation of the operator space. 

For that purpose, we first express the local operator A\ a by means of the creation and 
annihilation operators in the momentum representation as 

A= E alj{a\,_ a a k ,,- (J ){k\? l {k'\i)(i\k") . (Al) 

h h' h" 

Here (k\i) = (i\k)* = 1/VN exp(zfc • Ri). The operator a\ a 5{a\, _ a a k n _ a ) in Eq. (A-l) is the 
eigen state of Lq(z) = QL(z)Q, i.e., 

-e k » - S jfc //_ (T (2;)]|aJ. CT <5(4,_ (T a fe //_ CT )) . (A-2) 



Here T, ka (z) is the Fourier transform of T,ij a (z), which is defined by Y> kcr {z) = 
J2j ^jOa(z) exp(ifc ■ Rj). 

Substituting Eq. (A-l) into the approximate expression G^ a {z) = (A| CT | {z— Lq(z)^ 1 ^4] cr ) 
and making use of the relation (A-2), we reach the explicit expression of the screened cluster 
memory function in the RPT-0 as 

= E <#i><^i>< fc il*><^>(^ - (A-3) 

hh f h lf h-, h' h" 

=x(z-Lo-v c y 1 , (A-4) 

{X)k l k' l k'{c>kk'k"a = (4 1CT ' (5 ( a li-a' a fci'-^)l«L <5 ( a i'- CT a fe"-^)) » ( A-5 ) 
{Lo) klk ' lk '{cr'kk>k"a = (fcr + £fc + S fc(T (z) + efc' + Sfc'_o-(z) 

-Zk" ~'£'k"-a{z))&kik&k' 1 ki&k'{k"&a(j' > (A-6) 

{ v c) klk ' lk '{„i kk 'k"a = ^ [ - A; m{7 S; m(T (z)(fci|Z)(m|fe)(5 fc / fc /(5 fc // fc // 

Imgc 

- Aim-aSim-CT (z) ( ^ |Z) (m| k')5 klk 5 k » kl/ 

+\ lm -eE lm _ (7 (z)(k"\l)(rn\k")5k 1 kSk> 1 k' <W • (A-7) 

The screened memory function (A-3) depends on the choice of When we choose 

{Aijo- = 1} and adopt the Hartree-Fock approximation when computing the static average in 
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Eq. (A-5), we have 

n( c ) ST ( i \ k )( i \ k ')( k "\ i )x(ek*, ejs-a, ek>>-a){k\j)(k'\j)(j\k") 

u 0ija\ z ) — /_, ^ T\ o T~~ l — 7T > 

h k' h" Z ~ e kcr - ^ka(Z) - Z-W-o - ^k'-a(Z) + <i k » - a + l^ k n_ a {Z) 

X(e k , t w ,e v ,) = (1 - /(e fc ))(l - fMtfM + f {e k ) f {e k ,){l - /(e fc „)) • (A-9) 

Here e ka = e CT + e k is the Hartree-Fock energy, f(co) is the Fermi distribution function. 

A way to simplify in Eq. (A-8) the three-fold sum with respect to k might be to introduce 
an approximate T, ka (z) whose k dependence has been projected onto the Hartree-Fock energy 
e ka , as follows. 

Z a (e ka , z) = J \J ka> \ aK 1 . (A-10) 
J dk d(e ka - e k ' a ) 

After the approximation T, ka (z) « £ CT (efc CT , z), Eq. (A-8) is expressed as 

77(c) , v = /" dede'de' , pi j<T (e)p ij - <T (e')p j i- l r(e'')x(e, e', e") 
° 13 J z-e-t a (e,z)-e>-t_ a (e',z) + e" + t_ a (e»,z) ' 

Here Pij a (e) is the Hartree-Fock density of states defined by Pi J0 -(e) = X^(i|A:)c)(e ~~ e fco-)(fc| j)- 
On the other hand, in the case of { A^ Jcr = 0} (i.e., {Ay CT = 1}), it is not easy to obtain 

fc) 

directly a simplified expression of Eq. (A-3) because v c remains. However, (z) in Eq. 
(48) becomes the Coulomb interaction in this case. The second-order perturbation of the 
temperature Green function yields then an approximate expression. 

f dede'd^p^p^^p^y')^^) 
AZ) ~ J z-e-e' + 6" • (A ' 12) 

Here pf^ a (e) is the DOS for a cavity Green function for the Hamiltonian (47) with {Xij a = 1}. 

Therefore we obtain a simplified expression of the screened cluster memory function, which 
is an interpolation between Eq. (A-ll) for Ay CT = 1 and Eq. (A-12) for Ay CT = 0. 

{z j dede'de"p^(X, e)p^ a (X, e')p^(X, e") X (e, e>, e") 

^ i z-e-X a t a (e,z)-e> -X- a t_ a (e>,z) + e» + X- f7 £- a (e",z) ' 

This is Eq. (61). The prefactor Aij a , the densities of states for the cavity states p\f a {X, e), and 
a simplified self-energy S CT (e, z) are given by Eqs. (63), (65), and (67), respectively. 



77(c) 



77(c) 
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